Core percolation on complex networks 
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As a fundamental structural transition in complex networks, core percolation is related 
to a wide range of important problems. Yet, previous theoretical studies of core percolation 
have been focusing on the classical Erdos-Renyi random networks with Poisson degree distri- 
bution, which are quite unlike many real-world networks with scale-free or fat-tailed degree 
distributions. Here we show that core percolation can be analytically studied for complex 
networks with arbitrary degree distributions. We derive the condition for core percolation 
and find that purely scale-free networks have no core for any degree exponents. We show 
that for undirected networks if core percolation occurs then it is always continuous while for 
directed networks it becomes discontinuous when the in- and out-degree distributions are 
different. We also apply our theory to real-world directed networks and find, surprisingly, 
that they often have much larger core sizes as compared to random models. These find- 
ings would help us better understand the interesting interplay between the structural and 
dynamical properties of complex networks. 

Network science has emerged as a prominent field in complex system research, which provides us 
a novel perspective to better understand complexity [1^3j. In the last decade considerable advances 
about structural and dynamical properties of complex networks have been made[H-[6]. Among 
them, structural transitions in networks were extensively studied due to their big impacts on nu- 
merous dynamical processes on networks. Particularly interesting are the emergence of a giant 
connected component [7HIU] , fc-core percolation [TTH13] . fc-clique percolation[l51 [15], and explosive 
percolation |16l - fT8] . These structural transitions affect many properties of networks, e.g. robust- 
ness and resilience to breakdowns [91 [T9l [20]. cascading failure in interdependent networks [2TH23] . 
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epidemic and information spreading on socio-technical systems [3"1 [25], [26] . Recent work on network 
controllability reveals another interesting interplay between the structural and dynamical properties 
of complex networks |27H29| . It was found that the robustness of network controllability is closely 
related to the presence of the core in the network |27[ 130] . Actually, core percolation has also been 
related to many other interesting problems, including conductor-insulator transitions [HJ E2] and 
some classical combinatorial optimization problems, e.g. maximum matching [331435] and vertex 
cover [36H38]. 

The core of a undirected network is defined as a spanned subgraph which remains in the network 
after the following greedy leaf removal (GLR) procedure |32| [3"3~] : As long as the network has leaves, 
i.e. nodes of degree 1, choose an arbitrary leaf v\, and its neighbor V2, and remove them together 
with all the edges incident with t>2- Finally, we remove all isolated nodes. It can be proven that 
the resulting graph is independent of the order of removals [32 ]. Note that the core described above 
is fundamentally different from the fe-core of a network. The latter is defined to be the maximal 
subgraph having minimum node degree of at least k, which can be obtained by iteratively removing 
nodes of degree less than k. Apparently, the GLR procedure described above is more destructive 
than the node removal procedure used to obtain the 2-core (see Figflk). In studying the robustness 
of controllability for general directed networks, the GLR procedure has been extended to calculate 
the core of directed networks[27j. We first transform a directed network Q to its bipartite graph 
representation B by splitting each node v into two nodes v + (upper) and v~ (lower), and we connect 
Vi to in B if there is a link (v\ — > v?) in Q. The core of a directed network Q can then be 
defined as the core of its corresponding bipartite graph B obtained by applying GLR to B as if B 
is a unipartite undirected network. 

One can easily tell whether the core exists in two very special cases: (1) If a network has no 
cycles, i.e. a tree or a forest (a disjoint union of trees), then eventually all nodes will be removed, 
hence no core. For example, the Barabasi- Albert (BA) model with parameter m = 1 yields a tree 
network, hence no core exists. (2) If a network has no leaf nodes, e.g. regular graphs with all nodes 
having the same degree k > 1 or the networks generated by the BA model with m > 1, then the 
GLR procedure will not even be initiated, hence all the nodes belong to the core. 

Except those two special cases, no general rules have been proposed to predict the existence 
of the core for an arbitrarily complex network. Previous theoretical studies focused on undirected 
Erdos-Renyi (ER) random graph. It has been show that for mean degree c < e = 2.7182818. . ., 
the core is small (zero asymptotically), whereas for c > e the core covers a finite fraction of all 
the nodes [521 133] [39] . In other words, core percolation occurs at the critical point c* = e. More 
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interestingly, it has been suggested that in ER random graph core percolation coincides with the 
changes of the solution-space structure of the vertex cover problem |36l 1381 I40j . which is one of 
the basic NP-complete optimization problems [41]. Also, for c < e the typical running time of an 
algorithm for finding the minimum vertex cover is polynomial [32j [36], while for c > e, one needs 
typically exponential running time[l2]. Hence, core percolation also coincides with an "easy-hard 
transition" of the typical computational complexity [381 |4~U] . 

Despite the results on undirected ER random networks and the importance of understanding 
the intriguing interplay between core percolation and other problems, we lack a systematic study 
and a general theory of core percolation for both undirected and directed random networks with 
arbitrary degree distributions. 



I. ANALYTICAL FRAMEWORK 



We propose the following analytical framework to study core percolation on random networks 
with arbitrary degree distributions. We first categorize the nodes according to how they can be 
removed during the GLR procedure. We define the following categories: (1) a-removable: nodes 
that can become isolated (e.g. v\ and V2 in Figf!}}); (2) /3-removable: nodes that can become a 
neighbor of a leaf (e.g. v% and i>5 in FigjlJ)); (3) non- removable: nodes that cannot be removed 
and hence belong to the core (e.g. vq,vy and vs in Fig ^>). While the core is independent of the 
order the leaves are removed[32j, the specific way a node is removed may depend on this order, 
but it can be proven that no node can be both a-removable and /3-removable at the same time. 
Now we consider an uncorrelated random network with arbitrary degree distribution P(fc)[10l 143]. 
Assuming that in each removable category the removal of a random node can be made locally, 
we can determine the category of a node v in a network Q by the categories of its neighbors in 
Q \ v, i.e. the subgraph of Q with node v and all its edges removed, using the following rules: (1) 
a-removable: all neighbors are /3-removable; (2) /3-removable: at least one neighbor is a-removable; 
(3) non-removable: no neighbor is a-removable, and at least two neighbors are not /3-removable. 

Let a and /3 denote the probability that a random neighbor of a random node v in a network Q 
is a-removable and /3-removable in Q \ v, respectively. We can derive two self-consistent equations 
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about a and /3 

oo 

a = Y,Q( k )P k ~ 1=A (l-P), (1) 
fc=i 

oo 

1-/3 = Jg(*)(l-a) M = A(a) (2) 
fc=l 

where Q(/c) = kP{k)/c is the degree distribution for the node at a random end of a randomly 
chosen edge, c = X^fcLo kP(k) is the mean degree, and A(x) = J2T =0 Q(k + 1)(1 - x) fc . These two 
equations indicate that a satisfies x = A^A(x)^. It can be shown that a is the smallest fixpoint of 
A(A(x)), i.e. the smallest root of the function f(x) = A(A(x)) — x. 

The expected fraction of non-removable nodes, i.e. the normalized core size (n core = N core /N), 
can then be calculated: 

oo k / , v 

n cove = Y,P(k)Y,( )/3 fc " s (l-/3-a) s , (3) 

fc=0 s=2 

which can be simplied in terms of G(x) = YlV=o P{k)x k , i.e. the generating function of the degree 
distribution P(k). The final result is given by 

n COTe = G(l-a)-G(P)-c(l-l3-a)a. (4) 

For Erdos-Renyi random networks, G{x) = e _c ( 1_a: ) = A(l — x), Eqj4]can be further simplified as 
^core = (1 — /3 — a)(l — ca), confirming previous results [32"1 139], 

The normalized number of edges in the core (Z core = L core /N) can also be calculated in terms 
of a and /3. Consider a uniform random edge, which remains in the core if and only if both of its 
endpoints are non-removable without removing the edge. The probability of one endpoint being 
non-removable without removing the edge is 1 — a — (3, and for the two endpoints the probabilities 
are independent. Therefore, the expected normalized number of edges in the core is 

U e = |(l-«-/3) 2 . (5) 

with c/2 = L/N the normalized number of edges in the network. Clearly, both n core > and 
^core > if and only if 1 — f3 — a > 0. 

Now we consider directed networks Q with given in- and out-degree distributions, denoted by 
P~{k) and P + (k), respectively. Let c denote the mean degree of each partition in the bipartite 
graph representation B of the directed network Q, i.e. the mean in-degree (or out-degree) of Q. 
Define Q^(k) = kP^(k)/c, which is the degree distribution of the upper or lower end, respectively, 
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of a random edge in B. Define A (x) = YlkLo Q ft + — x ) k - Then the same argument as we 
used in the undirected case gives that 

a ± =A ± (l-^), (6) 
l-/3± = (7) 

and a 1 * 1 is the smallest fixpoint of (A^(x)) . Now we can calculate the size of the core for each 
partition in B as 

oo k / . \ 

ntore = £ £ ( J (/? T ) fc " S (l " /3 T " « T ) S (8) 

fc=0 s=2 ^ ' 

and we define the size of the core in the directed network Q as 

^■core = ( n core + ^corc)/^- (9) 

The normalized number of edges in the core can also be calculated 

Zcore = c(l-a + -/3 + )(l -aT-p-). (10) 

II. CONDITION FOR CORE PERCOLATION 

It is easy to see that the core in a undirected network with degree distribution P(k) is the very 
same as in a directed network with the same out- and in-degree distributions, i.e. P + {k) = P~ (k) = 
P(k). Therefore we can deal with directed network for generality. As n COTe is a continuous function 
of or* 1 , we focus on a^, which is the smallest root of the function / (x) = A^ (A^x)) — x. There 
are several interesting facts about the function /^(x). First of all, since A^(x) is a monotonically 
decreasing function for x 6 [0, 1] and ^^(O) = 1 is the maximum (see Figsj2jk3|), we have /±(0) > 
and f±(l) < (see Figj3j3,d). Consequently, the number of roots (with multiplicity) of f^(x) in 
[0,1] is odd, and numerical calculations suggest that this number is either 1 or 3 (see Figsj2j|3]). 
Secondly, if /^(xo) = then (^^(xo)) = 0, which means ^^(x) transforms the roots of /^(x) to 
the roots of / T (x). This also suggests that f (x) always has a trivial root or* 1 = ^(a^) = 1 — (5^ . 
(For undirected networks, /(x) always has a trivial root a = A(a) = 1 — f3.) Since ^^(x) is 
a monotonically decreasing function and is the smallest root of /^(x), A^(a^) = 1 — /3 T is 
therefore the largest root of /^(x). Hence 1 — — is the difference between the largest and 
the smallest roots of /^(x) (see Figj2|. Consequently, if /^(x) has only one root (which then must 
be the trivial root = A ± (a T ) = 1 — (3^), then 1 — (3^ — = 0. According to Eq{8j this 



6 



implies that there is no core. On the other hand, if multiple roots exist and they are different then 

1 — /3 — ot > 0, and the core will develop. 

We apply the above condition to the following random undirected networks with specific degree 

distributions[10j. (1) Erdos-Renyi (ER)[71[8] networks with Poisson degree distribution P(k) = 

e~ c c k /k\, A(x) = e~ cx and f(x) = exp(— ce~ cx ) — x. As shown in Figjj^,, the core percolation 

occurs at c = c* = e, which agrees with previous theoretical results E3 EH]- (2) Exponentially 

distributed graphs with P(k) = (1 — e _1 / K )e~ fc / K and mean degree c = e~ l / K /{\ — e _1//K ). We 

find that core percolation occurs at c = c* = 4. (3) Purely power-law distributed networks with 

P(k) = fc~ 7 /C(7) for k > 1, 7 > 2 and £(7) the Riemann £ function. We find that f{x) has no 

multiple roots and hence n core = for all 7 > 2. In other words, for purely scale-free (SF) networks, 

the core does not exist. (4) Power-law distributed networks with exponential degree cutoff, i.e. 

P(k) = (e- 1 /^) ^ or ^ — 1 with Li ra (x) the nth polylogarithm of x. We find that n core = for 

7 > 7c(«0, and the threshold value 7 c (k) approaches 1 as k increases. Hence, for SF networks 

with exponential degree cutoff the core still does not exist for all 7 > 1. (5) Asymptotically SF 

[-(l— f)] 1 ^ r(fe— £]) 

networks generated by the static model with P(k) = — — ^ r(fc+i) ' w ^ ere -^(s) is the 

gamma function and T(s, x) the upper incomplete gamma function [44H46] . In the large k limit, 
P(k) ~ k ^ 1+ f^ = k^" 1 where 7 = 1 + | > 2. For small k, P{k) deviates significantly from the 
power-law distribution [35] and there are much fewer small-degree nodes than the purely scale- free 
networks, which results in a drastically different core percolation behavior. 

Hereafter, we systematically study the net effect of adding more links (i.e. increasing mean 
degree c, yet without changing other parameters in P{k)) on core percolation. ER networks and 
the asymptotically SF networks generated by the static model naturally serve this purpose, since 
their mean-degree is an independent and explicit tuning parameter. 



III. NATURE OF CORE PERCOLATION 

We observed that if the mean degree c is small, then /^(x) has one root, but if c is large, f^{x) 
has three roots (see Figsj2j [3]). At the critical point c = c*, the number of roots jumps from 1 to 
3 by the appearance of one new root with multiplicity 2. (Note that f (x) cannot immediately 
intersect the x-axis at two new points, but it touches first.) This explains why the core percolation 
occurs at c = c*. 

According to the transformation from the roots of f^(x) to the roots of f^(x) through A^(x), 
for either f + (x) or f~(x) (depending on the details of P + (k) and P~(k)) its new root at c = c* 
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is smaller than its original root; and for either f~(x) or f + (x) the new root at c = c* is larger 
than the original root; or there is a degenerate case when this new root is the same as the original 
root for both f + {x) and f~(x). For example, for directed asymptotically SF networks generated 
by the static model with 7i n = 2.7, 7 ou t = 3.0, the new root (marked as green dot) of f + (x) at 
c = c* is smaller than the original root (green square) of f + {x) (see Fig^), and the new root 
(green square) of f~(x) at c = c* is larger than the original root (green circle) of f~(x) (see 
Fig|3ji). In other words, at the critical point, for either f + (x) or f~(x), its smallest two roots 
are the same, and for the other function (either f~(x) or / + (x)), its largest two roots are the 
same (see Figj3]:,d). While for directed networks with P + (k) = P~(k) = P(k), i.e. the degenerate 
case, we have f + (x) = f~(x) = /(x), and the new root of f(x) at c = c* has to be the same 
as the the original root of f(x), i.e. all three roots must be the same (see Figj3^,). Therefore at 
the critical point, unless in the degenerate case, a + together with j3~ (or a~ together with f3 + ) 
decrease discontinuously, which implies a discontinuous transition in the core size. To sum up, in 
the degenerate case that P + (k) = P~(k) = P(k) core percolation is continuous, but for general 
non-degenerate case P + (k) ^ P~{k), we have a discontinuous transition in both n CO re and ^ C ore- 
These results are clearly shown in Figj3]3,e. 

At the critical point c* , f^(x) touches the x-axis at its new root (see Figj3j:,d), hence we have 
either f + (a + ) = (/ + )'(a+) = (or - /T ) = - /T ) = 0), which enable us to calculate 

the core percolation threshold c*. In the degenerate case, if c < c* then f(a) = f'(a) = can 
be further simplified as A(a) = a and [A'(a)] 2 = 1. The results of c* for ER and SF networks 
generated by the static model are shown in Fig|4^. 

The discontinuity in n core and Z core at c*, denoted by A n and Aj respectively, can also be 
calculated 

A„=2( A n+An) (11) 

Ai = c*(l - /T'* - a~ *)(1 - - a+'*) (12) 

with A± ee G ± (l - a T '*) - G ± (/3^'*) - c* (1 - - a^'*)a^*. The results of A n for ER 
and SF networks generated by the static model are shown in Fig|4}3. We find that A n — > 
as 7i n —> 7out> consistent with the result obtained above that core percolation is continuous for 
undirected networks or directed networks with P + (k) = P~(k). We also find that A n increases as 
the differences between 7; n and 7 out increases. 

We can further show that in the general non-degenerate case, core percolation is actually a 
hybrid phase transition [T2], [131 HZ!) i- e - n core (or / C oro) has a jump at the critical point as at a 
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first-order phase transition but also has a critical singularity as at a continuous transition. The 
results are summarized here: in the critical regime e = c — c* — > + 

n core - A n ~ (c-cT (13) 
Ue-Ai ~ (c-c*) e (14) 

with the critical exponents 77 = 9 = g- Our calculations do not use any specific functional form 
of AF(x). Instead, we only assume that they are continuous functions of the mean degree c. 
Interestingly, in the degenerate or undirected case, one has a continuous phase transition (A n = 
Ai = 0) but with a completely different set of critical exponents: rf = 9' = 1 [32J. 



IV. NUMERICAL RESULTS 



We check our analytical results with extensive numerical calculations by performing the GLR 
procedure on finite discrete networks generated by the static model [4"4Tl46| . Figjj^i and [5)3 show 
n core and Z core (in symbols) for undirected ER networks and asymptotically SF networks with 
different degree exponents. For comparison, analytical results for infinite large networks are also 
shown (in lines). Clearly, core percolation is continuous in this case. This is fundamentally different 
from the k > 3-core percolation, which becomes discontinuous for ER networks and SF networks 
with 7 > 3PHH2!. 

FigjSJ: and[5]i show the results of n core and l COTe for directed networks. For directed networks with 
the same in- and out-degree distributions, e.g. directed ER networks or directed SF networks with 
7in = Tout generated by the static model, the core percolation is still continuous. But for directed 
networks with different in- and out-degree distributions, e.g. directed SF networks with 7i n 7^ 7 out 
generated by the static model, the core percolation looks discontinuous. The discontinuity in n core 
(or / core ) increases as the difference between 7i n and 7 ou t increases (see Fig|5kf). 



V. REAL NETWORKS 

We also apply our theory to real-world networks with known degree distributions. In Figj6] 
we demonstrate that in some cases our analytical results calculated from Eqsj4j [5] (or Eqs(9j |i~0| ) 
with degree distribution as the only input predict with surprising accuracy the core size of real 
networks. Yet, in other cases there is a noticeable difference between theory and reality, which 
suggests the presence of extra structure in the real-world networks that is not captured by the 
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degree distribution. In particular we find that almost all the directed real-world networks have 
larger core sizes than the theoretical predictions (see Figj6^i,b). In other words, those networks 
are "overcored" . While if we treat those networks as undirected ones, their core sizes deviate from 
our theory in a more complicated manner. The effects of higher order correlations (e.g. degree 
correlations |48j. clustering [39], loop structure [50] and modularity |51j) may play very important 
roles to explain the discrepancy between theory and reality. 



VI. CONCLUSION 



In sum, we analytically solve the core percolation problem in both undirected and directed 
random networks with arbitrary degree distributions. We show the condition for core percolation. 
We find it is continuous in undirected networks (if it occurs), while it becomes discontinuous or 
hybrid in directed networks unless the in- and out-degree distributions are the same. Within 
each case, the critical exponents associated with the critical singularity are universal for random 
networks with arbitrary degree distributions parameterized continuously in mean degree. But the 
two cases have totally different sets of critical exponents. These results vividly illustrate that core 
percolation is a fundamental structural transition in complex networks and its implication on other 
problems, e.g. conductor-insulator transitions, combinatorial optimization problems, and network 
controllability issue, deserves further exploration. The analytical framework presented here also 
raises a number of questions, answers to which would further improve our understanding of core 
percolation on complex real-world networks. For example, we focused on uncorrelated random 
networks and leave the systematic studies of the effects of higher order correlations as future work. 
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FIG. 1: The core of a small network, a, The core (highlighted in red) obtained after the greedy leaf 
removal procedure is fundamentally different from the 2-core (highlighted in green) obtained by iteratively 
removing nodes of degree less than 2. The 2-core contains the core, whereas the opposite is not true. Size of 
nodes are roughly proportional to the degree of nodes, b, Removal categories of nodes according to how they 
can be removed during the greedy leaf removal procedure. Red nodes are non-removable, i.e. they belong to 
the core. Green nodes are removable: nodes xi\ and V2 are a-removable; nodes v$ and v$ are /3-removable. 
White node is removable but it is neither a-removable nor /3-removable. Node v§ is /3-removable because 
t>4 will become a leaf node after removing node v\ (or V2) together with V3. 
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FIG. 2: Graphical solution of the self-consistent equations. a-d, For undirected networks, the 
function A(x) transforms the roots of f(x) to the roots of the same function f(x). The graphical solution of 
f(x) — A[A(x)j — x — is best illustrated by plotting the two curves A{x) vs. x (in red) and y vs. A(y) (in 
green) in the same coordinate system. The coordinates of the intersection point(s) of the two curves give the 
solution(s) of f(x) = 0. In a, b, and c, we show the graphical solutions for c <, =, and > c* , respectively, d, 
By drawing the two curves (A(x) vs. x) and (y vs. A(y)) at different mean degrees c, we get two surfaces. 
The intersection curve of the two surfaces yields the solutions of f(x) = at different c values. For c < c* , 
the intersection curve has one branch given by (a, 1 — j3, c) = (1 — /3, a, c). For c > c* , the intersection curve 
has three branches. The top and bottom branches are given by (a, 1 — /3, c) and (1 — j3, a, c), respectively, 
e-h, For directed networks, A ± (x) transforms the roots of f T (x) to the roots of f (x). The graphical 
solution of f ± (x) — A ± (A T (x)) — x — can be illustrated by plotting A + {x) vs. x (in red) and y vs. A~{y) 
(in green) in the same coordinate system. The x-coordinate (or y-coordinate) of the intersection point(s) of 
the two curves give the solution(s) of the equation f~(x) = (or f + (x) — 0, respectively). In e, f, and g, 
we show the graphical solutions for c <,=, and > c*, respectively, h, By drawing the two curves (A + (x) 
vs. x) and (y vs. A~(y)) at different mean degrees c, we get two surfaces. The intersection curve of the 
two surfaces yields the solutions of f (x) = at different c values. For c < c* , the intersection curve has 
one branch given by (or , 1 — /3 + , c) = (1 — a + , c). For c > c* , the intersection curve has three branches. 
The top and bottom branches are given by (a - , 1 — c) and (1 — (3~ , a + , c), respectively. 
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FIG. 3: Analytical solution of the core percolation, a-b, Undirected Erdos-Renyi (ER) random 
networks, a, a is the smallest root of the function f(x) = A(A(x)) — x, represented by red, green, and blue 
dots for c <, =, and > c* = e, respectively, b, a,f3,n corc and l COTe as functions of the mean degree c. c-e, 
Directed asymptotically scale-free (SF) random networks generated by the static model. Both the in-degree 
and out-degree distributions of the networks are scale-free with degree exponents 7i n = 2.7 and 7 ou t = 3.0. 
c, d, a is the smallest root of the function f ± (x) = A ± (^A T (x)^j — x, represented by red, green, and blue 
dots for c <,—, and > c* ~ 11.2, respectively, e, a*, {3 ± , n corc and Z corc as functions of the mean degree c. 
The jumps in a + and /3~ result in the jumps in n corc and Z C orc 7 hence the first-order core percolation occurs. 
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Yout 



FIG. 4: Threshold and discontinuity of core percolation, a, Analytical solution of the core percolation 
threshold c* calculated by solving / ± (a;) = f ± '(x) = for model networks. For ER networks, c* = e. For 
undirected asymptotically SF networks generated by the static model, c* — > oo as 7 — > 2, and and c* — > e as 
7 — > 00. b, The discontinuity A n in n corc at c = c* for model networks. For undirected or directed networks 
with P + (k) = P~(k), A n = 0. For directed network, A n increases as the difference between the in- and 
out-degree distributions (quantified by the difference between the degree exponents 7; n and 7 ou t) increases. 
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FIG. 5: Core percolation in random networks. Symbols are numerical results calculated from the GLR 
procedure on finite discrete networks constructed with the static model[44 with N = 10 5 . The numerical 
results are averaged over 20 realizations with error bars defined as s.e.m. Lines are analytical results for 
infinite large system (N — > oo) calculated from Eq|4] and [5] for undirected networks or Eq{9] and 10 for 
directed networks. Finite size effects arc more disccrnablc for 7 — > 2, which can be eliminated by imposing 
degree cutoff in constructing the SF networks |521 [55] . a-b, The normalized core size (n C ore = N cme /N) and 
the normalized number of edges in the core (7 C orc = L cmc /N) for undirected model networks: Erdos-Renyi 
(ER) and asymptotically scale-free (SF) with different values of 7. For both model networks, the core 
percolation is continuous, which is fundamentally different from the k > 3-core percolation, which becomes 
discontinuous for ER networks and SF networks with 7 > 3 [ 1 1 , 12J. c-d, n coro and Z coro for directed ER and 
asymptotically SF model networks. The core percolation is continuous if the out- and in-degree distributions 
are the same (P + (k) = P = {k)) while it becomes discontinuous if P + (k) ^ P = (k). c-d, For directed SF 
networks with fixed 7 ou t = 3.0, by tuning 7j n we see that the discontinuity in both n core and i cor e become 
larger as the difference between 7j n and 7 ou t increases. 
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FIG. 6: Normalized core size for real networks, compared with analytical predictions. All the 

real networks considered here are directed. For data sources and references, see Ref.[37] Supplementary 
Information Sec. VI. a-b, By applying the GLR procedure we yield n T ^ l e and l T c %f e - Using Eqj9]and Eq 10 



with out- and in-degree distributions (P + (k) and P~(k)) as the only inputs, we obtain n^^ ytlc and ^core ytlc - 
c-d By ignoring the direction of the edges, we can treat the original directed networks as undirected ones 
and apply the GLR procedure to get n T c %f e and ^or'- Similarly, we can calculate n™ a c lytlc and i^!^ by 
using Eqj4]and Eqj5]with the degree distribution P(k) as the only input. 
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